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Abstract 



We study the statistical mechanics of the self-gravitating gas at thermal equi- 
librium with two kinds of particles. We start from the partition function in the 
canonical ensemble which we express as a functional integral over the densities of 
■ the two kinds of particles for a large number of particles. The system is shown to 

possess an infinite volume limit when {Ni,N2, V) — > oo, keeping Ni/Vs and N2/V3 
O I fixed. The saddle point approximation becomes here exact for (A^i,A'^2)^) 00. 

' It provides a nonlinear differential equation on the densities of each kind of parti- 

es ■ cles. For the spherically symmetric case, we compute the densities as functions of 

two dimensionless physical parameters: m = '^'^/^^ and m = '^"j^'^^ (where G is 

- VST VST 

^ ■ Newton's constant, mi and m2 the masses of the two kinds of particles and T the 



temperature). According to the values of rji and r]2 the system can be either in a 
gaseous phase or in a highly condensed phase. The gaseous phase is stable for r]i 
and r]2 between the origin and their collapse values. We have thus generalized the 
well-known isothermal sphere for two kinds of particles. The gas is inhomogeneous 
and the mass M{R) inside a sphere of radius R scales with R as M{R) oc R"^ sug- 
gesting a fractal structure. The value of d depends in general on iji and r/2 except 
on the critical line for the canonical ensemble in the ??2) plane where it takes the 
universal value d ~ 1.6 for all values of N1/N2. The equation of state is computed. 
It is found to be locally a perfect gas equation of state. The thermodynamic func- 
tions (free energy, energy, entropy) are expressed and plotted as functions of r/i and 
r]2 ■ They exhibit a square root Riemann sheet with the branch points on the critical 
canonical line. The behaviour of the energy and the specific heat at the critical line 
is computed. This treatment is further generalized to the self- gravitating gas with 
n-types of particles. 
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1 Introduction 

The self-gravitating gases have remarkable physical properties due to the long range nature 
of the gravitational force. They are not homogeneous even at thermal equilibrium. This 
fundamental inhomogeneity character suggested that fractal structures can arise in a self- 
interacting gravitational gas0, ^. 

Self-gravitating gases are used to describe cold clouds in the interstellar medium as 
well as the large scale structure of galaxies. In both cases, self-gravitating gases provide 
scaling laws for the mass distribution with Haussdorf dimensions compatible with the 
observations ID, |^, ^. 

All particles have the same mass in the self-gravitating gases in thermal equilibrium 
considered till now both in the hydrostatic approach and in the statistical mechanics 
approach[Q. We study in the present paper the statistical mechanics of a non-relativistic 
gas with two kinds of particles with masses mi and m2- Such a system, besides its own 
physical interest, has obvious astrophysical motivations since cold clouds in the galaxy 
are formed typically by several kinds of particles. For example, hydrogen and helium. 

Let us begin by briefly recalling some results about the statistical mechanics of the 
self-gravitating gas with one kind of particles ||^ . It is a gas of non-relativistic particles 
of mass m interacting through Newtonian gravity. The gas is in a volume V and in a 
thermal bath at temperature T. In the usual thermodynamic limit {N, — oo and 
N/V is flxed) the gaseous phase is not stable and the system collapses in a very dense 
phase. In the dilute limit N, V oo and N/Vs fixed the system can exist in the gaseous 



phase. This is a dilute limit since the average density ^ goes as V 3 — > 0. The relevant 
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physical parameter of the system is r] = with Newton's constant G and the length 

L = V^. T] is the ratio of the characteristic gravitational energy ^ii^^Jl and the kinetic 
energy T of the gas. For 77 = the ideal gas is recovered. From t] = till a critical value 
?7o = 2.43450 . . . the gaseous phase is stable in the canonical ensemble. When t] reaches the 
value rjo the gas collapses in a very dense phase. The velocity of sound becomes imaginary 
at this point triggering unstabilities that lead to the collapse of the gas[Q. The saddle 
point approximation applies between 77 = and the point r]'^ = 2.517551 . . . (associated 
to the Jeans instability) where the determinant of small fluctuations is positive. At r]'" 
the determinant of small fluctuations vanishes and the saddle point approximation breaks 
down0]. Beyond rj'^ the gaseous phase is stable and the mean field approximation holds 
in the microcanonical ensemble. Solving the saddle point equation in the spherically 
symmetric case allows to obtain the particle density and the thermodynamic functions 
as functions of the physical parameter t]. As shown in ref. |Q the mean field approach 
is equivalent to the hydrostatic description provided the ideal gas equation of state 
is postulated in the latter approach. The mass distribution turns to exhibit a scaling 
behaviour function of i?[Hi. 

We consider in this paper the self-gravitating gas with two kinds of particles in the 
canonical ensemble. We recast the partition function as a functional integral over the 
densities of particles of each kind when the number of particles is large. The statistical 
weight for each configuration of densities turns to be the exponential of an 'effective 
action' which is proportional to the number of particles. Therefore, we can use the saddle 
point approximation in the thermodynamic limit to evaluate the partition function. The 
'effective action' turns to be the free energy as a functional of the particle densities. 

When the saddle point provides a minimum of the free energy, the density solution 
of the saddle point equation is the most probable. It is certainly exact for an infinite 
number of particles, since the minimized free energy exponentially dominates the partition 
function. That is, the mean field theory defined by the saddle point becomes exact for 
an infinite number of particles. The mean field approximation for the canonical ensemble 
ceases to be valid on a critical line in the (?7i,r/2) plane (see below). Beyond this critical 
line the saddle point is not a minimum of the free energy and it fails to reproduce the 
physics of the system. 

More precisely, we consider A^^i particles of mass mi and N2 particles of mass m2 
interacting through Newtonian gravity in a volume V and in a thermal equilibrium at 
temperature T. By analogy with the self-gravitating gas with one kind of particles, we 
consider the dilute thermodynamic limit: 

Ni, N2, — 00 keeping — 7- and — r fixed (1) 

The two relevant physical parameter are here 

GmlNi GmlN2 

where L = V^ for a cubic geometry. Notice that eq.(|lD implies that the ratio N1/N2 stays 
fixed for iVi, iVg, ^ 00. 

The self- gravitating gas with two kinds of particle behaves as a perfect gas in the 
extremely diluted limit r/i ^ and ri2 0. When rji and/or ri2 grow, the gas becomes 
denser till it collapses into a very dense phase when rji and 172 reach the collapse line for 
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the canonical ensemble in the (?7i,//2) plane. By analogy with the gas with one kind of 
particles |^ we expect the collapse line to be very close and below the critical line in the 
{fli-iV^} plane. The gaseous phase keeps stable in the microcanonical ensemble beyond 
this critical line. 

We find here that the saddle point equations are two coupled non-linear differential 
equations for the densities of the two kinds of particles: pi (x) , p2 (x) . We succeed to reduce 
these equations to a single non-linear differential equation. We solve it in the spherically 
symmetric case. We express the densities as functions of the physical parameters r/i and 
?72- We thus find the isothermal sphere with two types of particles. 

We compute the mass inside a sphere of radius R centered at the origin and show that 
it scales with a Haussdorf dimension d. d decreases with rjf and r^^ from the value c? = 3 
for the ideal homogeneous gas till d ~ 1.6 in the canonical critical line. The Haussdorf 
dimension keeps decreasing beyond the canonical critical line in the stable phase of the 
microcanonical ensemble, d at the canonical critical line turns to be independent of the 
ratio ri^/ri2 and coinciding within the numerical precision with d ^ 1.6 for the canonical 
critical point of the gas with one kind of particles |Q. This indicates that the Haussdorf 
dimension at the canonical critical line is an universal value, independent of the gas 
composition. 

The dependence of the critical values of the parameters rji and 772 with the number of 
particles ratio is computed. The thermodynamic functions (free energy, energy, entropy, 
local pressure and pressure contrast) are expressed as functions of the physical parameter 
r]i and r]2- The pressure contrast [ratio between the pressure at the origin and the pressure 
at the boundary] turns to be lower for this mixture of particles than for the gas with one 
kind of particles. 

We compute the pressure at a point r of the gas and show that it locally obeys the 
equation of state of an ideal gas, 

P(r) = -^pi(r) + -^p2(r). (3) 

Since the gas is inhomogeneous, the pressure acting on any non-infinitesimal volume of 
the gas does not obey the ideal equation of state. (This was already the case for a gas 
with all particles of equal mass[§|). We plot in figs. the pressure of the gas at the 
surface. 

The mean field equations have a straightforward hydrostatic interpretation. We show 
that the mean field equations derived from the partition function are equivalent to the 
hydrostatic equilibrium equations provided that the ideal equation of state is postulated 
in the latter approach. We stress that we give here a microscopic derivation of the 
equation of state (|^) from the partition function. 

We then consider a self-gravitating gas formed by n kinds of particles with different 
masses. The mean field equations are derived and shown to reduce to a single non-linear 
differential equation. 

This paper is organized as follows. In section II we present the statistical mechanics of 
the self-gravitating gas with two kinds of particles in the canonical ensemble, in sec. HI 
we present the main thermodynamic magnitudes, the equation of state and the scaling 
behaviour of the particle distribution for spherical symmetry. In sec. IV we present the 
generalization for n kinds of particles. The Appendixes contain relevant mathematical 
developments and the hydrostatic approach to a self-gravitating gas with two kinds of 
particles. 
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2 Statistical mechanics and mean field theory for the 
self-gravitating gas 

We present here the partition function for the self-gravitating gas with A^i particles of 
mass mi and N2 particles of mass m2 inside a finite volume V and derive the mean field 
approach to it. 



2.1 The canonical ensemble 

We study the statistical mechanics of the self-gravitating gas with two kinds of particles 
in the canonical ensemble. The Hamiltonian of the system is 

^ G ^ G ^ G mim2 

l<i<j<Ni l^l.i ~ ^Ijl l<i<j<Af2 1^2'' ~ l<j<Afi,l<i<A^2 1^^'' ~ ^^jl 

Here, pi_i and qi^i are the momenta and the coordinates of the particles of mass mi. p2,i 
and q2,i are the momenta and the coordinates of the particles of mass m2. Therefore, the 
classical partition function of the gas is 

It is convenient to introduce the dimensionless coordinates variables ri^ = L (\\ \ and 
i'2,i = q2,i- The momenta integrals are computed straightforwardly. Hence the partition 
function becomes the product of the partition function of perfect gases with masses m\ 
and m2 times the coordinate integral Zint- 



3N2/2 



Zi 



int ; 



where 



„ „ l=Ni k=N2 

Zint =/.../ n '^^''1.1 n exp (?7i Mil + V2 U22 + y/ViV2 U12) 

'' 1=1 k=l 



and 



1^1 11 

^11 = 1^ 2^ 1 r ' "22 = ^ 2^ -. 

VN1N2 l<i<N,,l<j<N2 l^l.i ~ ^2,j| 



2.2 Mean field theory 

We approximate the function Zint ioT a large number of particles (A'^i ^ 1 and A^2 ^ 1) 
generalizing the approach of ref. [H, |^ for two kinds of particles. The function Z^t is 
expressed as a functional integral over all configurations with particle densities pi(x) and 
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P2(x) [pi(x) stands for the density of the particles of mass mi and P2(x) stands for the 
density of the particles of mass 7712]. 



Zint = J Dpi{.) Dp2{.) dbi db2 



with, 



F-Fa 



~ y J \x-y\ P^^^' P^^y^ ~ J |x-y| ^^^^^ ^'^^^ 

I [ d^xd'^y 

-JN1N2 y^rnrn / n r /^iW ^2(y) + 

|x-y| 

+A^i J d^x pi(x) Inpi(x) + y" d^x p2(x) lnp2(x) + 
+iNi bi[l- f d^x pi(x)] + iN2 62(1 - / d=^x p2(x)] . 



r]2 f d^x d'^y 



Fq = -NiTln 



eV / rriiT 

iVi 



iVoT In 



eV /m2TY 



(4) 



(5) 

bi and 62 are Lagrange multiplier enforcing the normalization of the densities: 

y'd3xpi(x) = i , y'd3xp2(x) = i . (6) 

F stands for the free energy of the gas for the pair of densities (pi,p2), while 



is the free energy of the perfect gas with masses rui and 7712- One recognizes the gravita- 
tional energy in the first two lines of eq.(BD, while the third line contains the entropy. 



Since the free energy becomes large in the thermodynamic limit {Ni,N2 ^ 1), the 
functional integral Zint @ is dominated by the minima of Extremizing the free 

energy with respect to the pair of densities (pi,P2) yields the saddle point equations 



f d'^y f 

Inpi(x) = ai + T]i :pi{y)+fJ'V2 

J V — X J 



lnp2(x) 



y-x| 
dV 



|y-x| 



-P2(y) 



(7) 



a2 + - Vi 

P J |y-x 



■pi(y) +V2 I I ^ ^ | P2(y) • 



|y-x| 



where we used eq.(|]). These equations define the mean field approach. We set ai = 
— 1 + ibi and 02 = —1+262 and we denote by p the ratio of masses of the two kinds of 
particles, 



mi 

777,2 



We set, 



pi(x) = exp[<l>i(x)] , p2(x) = exp[$2(x)] . 
Eqs.(0) give for the gravitational potential 



T 

nil 



[$i(x)-ai] 



T 

7Tl2 ■ 



[$2(x) - 02] 



(9) 
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created by the matter densities pi(.) and p2(-) at the point x. Using eqs.(|D-@, we see 
that these densities obey the Bohzmann laws 

pi(x) = e"^ e-^^W , p2(x) = e'^^ e'^^W , (10) 

containing the energy of a particle in this mean field gravitational potential, as it must 
be. and e"^ play the role of normalization constants. 

Applying the Laplace operator to the saddle point equations (|^) we find the differential 
equations 

A$2(x) +47r ^ e*iW +47rr/2 e*^W = 0. (11) 
^^ 

These equations are scale covariant. If ($i, $2) is a pair of solutions of eqs.([TT|), then the 
pair ($ia,$2a) defined by 

$ia(x) = <l>i(Ax) + lnA2 , $2a(x) = $2(Ax)+lnA' (12) 

is also a solution of eq.(|lT]). This property is due to the scale behaviour of Newton's 
potential. Using eq.(^, we reduce eqs.(0) to a single equation 

A$i(x) + 47r 7]i e*i('') + 47rp r]2 e"'"^ = . (13) 

Using eqs.(|) and ([T0|), eq.(^3|) becomes in dimensionless coordinates, 

Af/(x) = [mi iVi pi(x) + m2 P2(x)] . (14) 

We show in Appendix B that this equation is the condition of hydrostatic equilibrium for 
a two-component fiuid once the ideal gas equation of state is postulated. 

Therefore, the mean field approximation is equivalent to the hydrostatic description 
in the gaseous phase provided the ideal gas equation of state is assumed in the latter 
approach. 

Notice that local equations of state other than the ideal gas are often assumed in the 
context of self-gravitating fiuids@. As stressed in ref.0], one needs long range forces 
other than gravitational in order to obtain a non-ideal local equation of state in thermal 
equilibrium. 



3 Spherically symmetric case 

For spherically symmetric configurations the mean field equations become ordinary non- 
linear differential equations. We express here the various thermodynamic quantities in 
terms of the solution of a single ordinary differential equation. Such equation reduces to 
the well-known isothermal sphere if all particles have identical mass. 

3.1 Reduction of the equation of saddle point 

We consider here the spherically symmetric case where the mean field equations ([TT|) take 
the form 

1^ + 1^ + 4vr r^i e*^(«) + 4vr p r^2 e*^(^) = (15) 

d2<l>2 2 d$2 , . Til , . n 



where we work in a unit sphere. Therefore the radial variable runs in the interval < 

R < Rmaxi Rmax = ^ ■ Usiug the scale covariance of the mean field eqs. ([TT|) by the 
transformation ([l2|) , we can set 



Rmax ) \ S"/]! 

Rmax / \ 3?72 



MR) = X2(A^)+ln|^ . (16) 



We use the new parameters r^f = rji/ R^ax and r/f^ = rj2/ Rmax in the spherically symmetric 
case. Hence the mean field eqs. (^) are transformed into a reduced system of the form, 

Xi(A) + ^x'i(A) + e>^^W+/ie>^^W = 

X2(A) + f x;(A) + -e>^^W + e>^^W = 0. (17) 

Let us find the boundary conditions of these equations. In order to have a regular solution 
at origin we impose 

x;(0) = , X2(0) = 0. 
The system (|T^) is invariant under the transformation 

A ^ A , Xi ^ Xi - 2a z = l,2. 

Hence, we can choose 

Xi(0) = 0. (18) 

without loosing generality. As we see below, the remaining boundary condition on X2(0) 
is not independent from rji and ri2. 

The densities of the two kinds of particles pi and p2 are to be normalized according 
to eq.d^). We obtain using eqs.(||) and ([iBD , 

r^f = \ l^dx e^^(") , ri2=\ t e"^'^"^ • (19) 
A JO A JO 

Using the reduced mean field equations (|l3), it is straightforward to show from eq.(|19D 
that 

< + = -Ax'i(A) (20) 

+ = -AX2(A). 

Hence 

x;(A) = - xUA) . (21) 
Recalling that U is the gravitational potential (||) we see using eqs. ([l6D , that 

— x'i(A) = — x;(A) 

mi 7712 

is the gravitational field at the boundary of the sphere (i? = Rmax) consistently with 
eq. (|2lD . From eq.(^) we introduce a new parameter 

C^X2(A)--Xi(A) (22) 
8 



independent of A and function only of the physical parameters rj^ and 77^. Notice the 
boundary condition 

X2(0) = C . (23) 



The reduced mean field equations ([T7|) become a single equation with the parameter C 
as coefficient 

+ Y x'i(A) + e^^W + /i e^^^^^^ = (24) 
A 

and the boundary conditions Xi(0) = 0, x'i(O) = 0. 

Using eqs-d^) and (0) we can express the densities of the two kinds of particle in 
terms of the solution of eq.( 



\ 2 \ 2 o 

p^(R) = A_ e^'(^n^) , p^{R) = 2_^eX^^) . (25) 

I2 



3rii ' 3ri2 



Both Xi('^) s-iid X2(A) are functions of C [see (|2^)-(p^]. Inserting Xi{^) and X2{x) in 
eg ■ (p!9D , we see that r^f and r^^ become functions of A and C. Then expressing A and C 
as functions of r]^ and r^^, both densities of particles in eq.(p5D become functions of the 
radial variable R and of the physical parameters r]^ and 7]2- 
We find asymptotically from eq.(^) 

Xi(A) '=°°-2/i logA + 0(l) 

for /i > 1. Notice that this asymptotic behaviour may apply for non- physical values of A 
where the gas is actually collapsed. 

We now compute the thermodynamic quantities as functions of the physical parameters 
rji and ri2- 

3.2 Free energy 

Let us start by computing the free energy. Using eqs. (^-(||) we find 
T 2 



ai + 1 d=^x $i(x) e*^(")) + T (""^ + / '^2(x) e*^w) . (26) 

We express now the Lagrange multipliers ai and 02 as functions of the physical parameters 
77^ and T]2- In the spherically symmetric case the integration over the angles in eqs.(^ 
yields 

$i(i?) = a,+ 4vrr/i (^i J^"" dR' R'' e*^(«') + J^^^^ dR' R' e*^(«') j 

+4vr /. V2 di?' i?'^ e-^^^^') + di?' i?' e*^(^')) 

$2(i?) = a^ + An^dr dR' it^ e*^(^') + /'''""^ di?' i?' e*^^^')') 
fi \RJo JR J 

+Anr^2 (-^ Z^"" di?' e*^(«') + /J''""^ di?' i?' e*^(^')^ . (27) 
We introduce the densities for the two kinds of particles at the boundary (i? = Rmax) 

^ ^ gXi(A) ^ ^ eMi?™..) = ^ (28) 

9 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Figure 1: Free energy versus r/^ for ^ = 4, ^ = | and therefore \ 
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where we used eqs.(§) and ([TBI). Notice that A and C are functions of 1]^ and 7]^ as 
explained by the end of sec. 3.1. Hence, /i and /2 are functions of r]^ and r^^. 

We find for R = Rmax the following expressions for the Lagrange multipliers using the 
normalization of the densities flBD 



ai = In/i - r^f - /i?7^ 



02 = In /a - - r^f - r/^ 



(29) 



Inserting these expressions (E^) into the free energy (EBf), we find 



2 



AT, 



ln/2--r7f-r7f 



j d^x $i(x) e*^W + ^ /d'x <l>2(x) e*^W . 



(30) 



We compute the integrals in the second line in appendix A and we find 



1 



Ari[ln/i - r/f - /i 7^2^ + 3(1 - h)] + iV2[ln/2 - - < - + 3(1 - /2 



The free energy as well as the other physical quantities are functions of rjf and ri2- The 
parameters r^f and are linked by the relation, 



V2 



ml N- 



2 R 

'11 



mf Ni 



Hence, for fixed ratios and the physical quantities depend only on 7]^ or on r^^. 
In that case it is simpler to see the physical quantities as functions of r/f^ or r/f on a two 
dimensional plot than to watch the three dimensional surfaces for the physical quantities 
as functions of and rj^. 

We plot in fig.0 the free energy as a function of r/^ for fixed ^ = y. In the limit 
where the particles of mass mi dominate (A^i 3> N2) the free energy becomes 



T 



Ari[ln/i -r^f + 3(1-/1)] . 



We recognize here the free energy of a self-gravitating gas with A'^i particles of mass mi 
(see ref.H). 



3.3 Energy 

We compute here the gravitational energy of the gas. The density of gravitational energy 



IS 



1 



ep X 



/ miA^ipi( 



2^ ^ + ^^^2^)c/(x) 



where U is the gravitational potential (pi). Hence, 



ep X 



e*^W [$i(x) - ai] - ^ e*^W [$2(x) - a,] 



2V 



2V 



(31) 
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Using the expressions for the Lagrange multiphers (pQ]), we obtain for the gravitational 
energy density in the spherically symmetric case 



ep X 



NiT 

w 



2V 



Pi{R) 

P2{R) 



In 



' pl{Rn 



In 



Pi{R) 

P2{Rm 



P2{R) 



P 



-Vi - V2 



Integrating the energy density (^Tj), with the help of eqs. and we obtain 



NiT 
2 

NiT 



In /i - - r]i 



+ 



1 



P 



J d^x $i(x) e*iW - ^ jd^x <l>2(x) e*^W 



(32) 



(33) 



2 7 ^ ' 2 

The integrals in the second line are computed in appendix A yielding for the gravitational 
energy 



Ep = 3T[N,{f,-l) + N2{f2-l)] . 



(34) 



3.4 Entropy 

Using the gravitational energy ( |5^ and the free energy ( pOp we obtain for the entropy 

5 = 5o + iVi [6 (A - 1) - In A + r/f + /i 7^2^] + iV2 [6 (/2 - 1 ) - In /2 + i r]i + V2] 

P 



where 5*0 is the entropy of the perfect gas. 



So = N,\ In 



V /miTY 

iVi 



N,. In 



V /m2T\2 

iV2 



+ 



We plot in Fig. || the entropy S' as a function of r^l*^ for ^ = y . 



3.5 Local pressure 

Since the system is non-homogeneous the local pressure is not uniform. The density of 
gravitational force is 



F(r) 



mi Ni pi(r) + m2 N2 P2(r) 
V 



^radff/ (r)] 



where U(r) is the gravitational potential (H). 

We obtain for the force at the point r using eqs.(|I) and (|^), 

F(r) = ^ grad(e*i('-)) + ^ grad(e*^('-)) 

The link between the density of force and the pressure is 

F(r) = grad[P(r)] . 
Using eq.(|D, the local pressure is given by, 

P(r) = ^Pi(r) + ^P2(r) 



V 



V 



(35) 
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This is the local equation of state for the self-gravitating gas with two kinds of particles. 
We see that it locally coincides with the equation of state of a perfect gas. Since the gas is 
inhomogeneous, the pressure acting on any finite volume will not obey the ideal equation 
of state. 

For a point r at the boundary, using eqs.(|), p8| ) and (^) yields the external pressure 
function of t]^ and 77^ as, 

P = -^/i + -^/2 , (36) 

where /i and /2 are defined by eq.(^). This is the equation of state of the gas as a whole 
that we plot in fig.^. 

Combining eqs.(0) and ( PB| ) yields the virial theorem 

PV Ep 

= Ni + N2 + — . 

T 3T 



3.6 Physical behaviour of the system 

This self-gravitating system formed by two kinds of particles can be in two phases: gaseous 
and condensed. The first one corresponds to 7]i and 7]2 between the origin and their 
collapse values. In the gaseous phase the free energy has a minimum for the pair of 
densities (pi,p2), solutions of the saddle point equations (0). This pair of densities is the 
most probable distribution which become absolutely certain in the thermodynamic limit. 
All thermodynamic quantities follows from this pair of densities (pi,p2)- 

In the condensed phase, {pi, P2) from mean field does not describe the particle dis- 
tribution and the mean field approach fails to describe the condensed phase. It may be 
studied by Monte-Carlo methods as in ref.[|]. 

When the physical parameters t]^ = = we retrieve the perfect gas. When 7]^ and 
r]2 increase, the gas becomes denser at the center of the sphere {R = 0) and less dense at 
the boundary {R = Rmax) because of gravitational attraction [see Figs. |^^. This effect 
is more acute for the heavier particles showing that more massive particles diffuse to the 
denser regions. 

The equation of state is depicted in fig.^. In the ideal gas limit, rji = r]2 = 0, PV = 
{Ni + N2)T. In the case 77^ = (gas of particles of mass 7/12) and in the case r/f^ = (gas 
of particles of mass mi) we recover the equation of state of the self-gravitating gas with 
one kind of particles 0]. 

We call critical values ?7f" and r^^*" the points where PV/[{Ni + N2)T] exhibits a 
vertical slope, r^f-^ and 7)2'" define a critical line in the {t]^, 77^) plane. Namely, for each 
value of N1/N2 we have a different pair of critical points t]^'" and 77^*". We plot the critical 
line in the (7/f , 7/|*^) plane in fig.^. 

The surface pressure has a rim on the canonical critical line [see fig. ^. The projection 
of this rim on the {rj^, r/^) plane yields the critical line plotted in fig.^. For a fixed N1/N2 
we get a section of the equation of state surface depicted in fig.|[ This section turns to 
have a form analogous to the equation of state for the self-gravitating gas with one kind 
of particles @] . 

By analogy with the gas with one kind of particles, we expect that the gas collapses 
in the condensed phase for values of r]^ and 77^ slightly below r/f*" and 77^'" where the 
saddle point approximation breaks down. 
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PV/[(N1 +N2)T] vs. eta1 R and eta2R 
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0.5 1 1.5 2 2.5 etal R 3 

Figure 4: The critical line in the [r]^, r^^) plane. The mean field approximation is valid 
in the region below the critical line 



The physical quantities exhibit a square root Riemann sheet structure as functions of 
7]^- and ?7^. The branch points are on the critical line for the canonical ensemble. The 
lower branch [see Figs. ^ and |^ for the equation of state and figs. ^[]^ for the densities] 
describe a phase absent in the canonical ensemble. Such phase is realized and is stable in 
the microcanonical ensemble as it was the case for the gas with one kind of particles^. 

We plot in Fig.^the value of ^^f +/i at the critical line as a function of the number of 
particles N1/N2. This quantity is proportional to the total mass of the gas mi Ni+m2 N2. 
This critical parameter T^f*" + fi rj2^ interpolates between the two limiting cases (iVi ^ N2 
and A^^i ^ A^2) where one kind of particles dominate over the others. When A^^i ^ A^2 the 
particles of mass mi dominate and t]^^ ^ ■ this limiting case, rj^^ + ^ r^l^*^ — 
^RC _ 2.518 . . . which is the critical value for the self gravitating gas with one kind of 
particles associated to the Jeans instability |P, When N2 ^ A''i the particles of mass 
m2 dominate and 77^*" 3> rji'" . In this limiting case, r/f" + /i r/^*" ^r]^'-'' = 10.04 . . . for 
the mass ratio /i = 4 corresponding to a mixture of hydrogen and helium. 

From eq.(^) we see that the partial pressures of the particles with mass mi and m2 
are given by, 

Pi(r) = ^e*^W , P2(r) = ^e*^W. (37) 

The pressure contrast is defined by the ratio of the pressure at the center and the pressure 
at the boundary: P{0)/P{R^ax) §• We find from eqs.(0), (^, and 

P(0) l+zi^e^ 

a = — 



P{Rma.) eXiW+^2eC+Xi(A)/M " 

We extend this definition to each kind of particles and say that the partial contrasts are 
given by Pi{Q) / Pi{Rmax) for the particles of mass mi and -P2(0)/P2(-Rmax) for the particles 
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Figure 6: Critical value of the physical parameter rji + jj, versus In for — = 4. 



of mass m2. We find from eqs. ([T8|) , (0), (|25| ) and 

a,^-^^ = e--W and ^ = e--(^)/- . (38) 

We plot the contrast and the partial contrasts in fig. ^. We see that the contrast a takes 
here lower values than for a gas with one kind of particles. On the contrary, ai, the 
partial contrast for the heavier particles, takes higher values than the contrast for a 
gas with one kind of particles. This is due to the fact that the overdensity of particles in 
the center is more acute for the heavier ones as noticed above. 

The particle density and pressure has its maximum at the origin and its minimum at 
the boundary, as expected. However, their ratio (contrast) is much larger for the heavier 
particles than for the lighter ones. They are related through [see eq.( 



where /i = mi/m2. Since 0:2 > 1, when mi > m2, we can get ai ^ ^2 [see fig.|||. 
In particular, ai may be much larger than the contrast in the gas with one kind of 
particles H, In summary, this shows that the heavier particles diffuse to the denser 
regions. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



Figure 7: Density of particles of mass mi [pi{Rmax)] and density of particles of mass 
m2 [p2{Rmax)] at the boundary versus where the mass ratio is ^ = 4, the number of 

particles ratio is ^ = | and then \ = ^. 



3.7 Scaling law 

We compute here the mass M{R) inside a sphere of radius R {0 < R < Rmax)- 

Using Gauss's theorem and recalling that U is the gravitational potential (|^) we find 
that 

M{R) = -^R' 

Using eqs.(0) we obtain 

'11 ^ -'^max ' ^ ^raax ' 

As for the self-gravitating gas where all particles have the same mass, the mass M[R) for 
the self-gravitating gas with two kinds of particles follows approximately the scaling law 

M{R) ^CR'^ . 

This indicates a fractal distribution with Haussdorf dimension d. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Figure 8: Density of particles of mass mi [pi(0)] and density of particles of mass 
"^2(0) [p2(0)] at the origin versus ry^ where the mass ratio is ^ = 4, the number of 

particles ratio is ^ = | and then ^ = ^. 

d decreases with 77^ and 77^^ from the value d = 3 for the ideal homogeneous gas 
{r]i = = 0) till d ^ 1.6 in the canonical critical line. The Haussdorf dimension keeps 
decreasing beyond the canonical critical line in the stable phase of the microcanonical 
ensemble. 

We plot in fig.^ the mass M{R) for M{R) greater than 10% of the total mass of 
the gas for several values of r/f choosing the ratio //f /r]^ to be 0.534 . . .. We exclude the 
region M{R) < 0.1 where the mass distribution is almost uniform. This local uniformity 
is simply explained by the absence of gravitational forces at the origin R = due to the 
spherical symmetry. 
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Figure 9: Contrast P (0) / P (Rmax) , partial contrast Pi (0) / Pi (Rmax) for the particles of 
mass nil and partial contrast P2{0) / P2{Rmax) for the particles of mass m2 versus rj^ 
where the mass ratio is ^ = 4, the number of particles ratio is ^ = | and then 

= ^. We have at the critical point Pi (0) / Pi (Rmax) ^ 81. > P{0)/P{Rmax) ^ 9.9 > 

P2iO)/P2iRmax) ^ 3.0. 



TABLE 1. The fractal dimension d and the proportionality coefficient C from a fit 

Ni 
N2 



according to M{R) ^ C i?*^ for ^ = 4, ^ = 0.0334 and then % = 0.534 and different 



m2 



values of r/f and 



R 



As shown in fig. |TT|, the Haussdorf dimension at the canonical critical line turns to 
be independent of the ratio f]i/f]2- d coincides within the numerical precision with the 
Haussdorf dimension c? ^ 1.6 at the canonical critical point of the gas with one kind of 
particles 0]. This indicates that the Haussdorf dimension at the canonical critical line is 
an universal value, independent of the gas composition. 
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Figure 10: In M(R) versus In ^ for ^ = 4, ^ = 0.0334. . . and then % = 0.534. . . 

for different values of ?7f : rjf- = 0.01, r]^ = 0.03, rj^ = 0.06 and the canonical critical 
point r/f^ = 0.8002 . . . 







d 


C 


0.03627 


0.0187 


1.62 


1.19 


0.6682 


1.688 


1.61 


1.04 


1.096 


1.126 


1.65 


1.03 


1.415 


0.7978 


1.62 


1.04 



TABLE 2. The fractal dimension d and the proportionality coefficient C from a fit 
according to M{R) ^ C on the critical canonical line for different values of r^f-^ and 



3.8 Critical behaviour of the thermodynamic functions 

According to the behaviour of the free energy near the critical line (see Fig. 0), we find 
that the first derivatives of the free energy (energy, pressure) are continuous, while the 
second derivatives (specific heats, compressibility) are discontinuous. Using eq.(plD and 
the form of the functions /i and /2 near the critical line (see Fig. we find that the 
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Fieure 11: In M(R) versus In „ for — = 4 on the critical canonical line for different 
values of r/f and r]^: r/f^ = 0.03627. . . , r/^^ = 2.486 ... ; r/f ^ = 0.6682 . . . , r]^^ = 
1.688 ... ; r/f^ = 1.096 . . . , r]^^ = 1.126 ... and <^ = 1.415 . . . , t]^^ = 0.7978 ... The 
Haussdorf dimension d turns to be independent of the composition of the gas taking the 
value d ^ 1.6. It coincides within the numerical accuracy with the Haussdorf dimension 
for the gas with one kind of particles at the canonical critical point [Q. 



energy has two branches E_^. and E_ which behave near the critical point as, 

E± = liN, + N^)T ± DT^^i^^ (39) 

where D is a positive constant. Deriving the energy (|39D with respect to T and recalling 
eq.(^) for 772, we obtain the two branches of the specific heat at constant volume 

S DT 

Here, Ej^ and stand for the gaseous stable phase in the canonical ensemble and 
and Cy- correspond to a phase only realized in the microcanonical ensemble @]. We 
see that = E^ at the critical point and hence E is continuous at criticality while 
exhibits there an infinite discontinuity. (As is clear, negative values of cannot be 
realized in the canonical ensemble^). 
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4 A self-gravitating gas with n kinds of particles 

The generahzation of the treatment given in previous sections to a self-gravitating gas 
with n kinds of particles is straightforward. We give below the mean field equations and 
the more relevant results. 

The relevant parameters of the gas are now 

Gm^Ni ^ . 
Vi = j^rj, , 1 < « < ^, 

where Ni is the number of particles of mass m^. We assume that Ni/L stay fixed while 
L oo, Ni — s> 00,1 < i < n. Therefore, the ratios Ni/Nj also stay fixed for L 
00, A^j ^ 00, 1 < i, j < n. 

The coupled mean-field integral equations for the densities of particles pj(x) (1 < z < 
n) take the form 

lnp,(x) =ai + mi J2— [ Pi(y) , 1 < « < ^ , (40) 

tr^ m,- J y - X 



j=i 



where Lagrange multipliers ai, a2, a„ enforce the normalization of the densities. 

y d^xpi(x) = 1 , l<i<n. (41) 
Eqs.(l40|) give for the gravitational potential 



f/(x) = [$i(x) - ai] = -— [$2(x) - as] = ... = -— [$„(x) - a„] . (42) 
nil ^2 fTT'Ti 

Setting, 

a(x) =exp[<l>i(x)] , l<i<n (43) 



and applying to eqs. (|40|) Laplace operator, we find the partial differential equations 

A^i{x) + AnmiY ^ e*^^''^ = , l<i,j <n. (44) 
Using eq.(^2p, we reduce eqs.(^) to a single equation 

A$i(x) +47rmiJ2 — e"^"^"^ e^*^^'^) = . 



Eqs. ( PI) are scale covariant. If $1, $2, •••5 are solutions of eqs. (|^), then $1^, ^2X, ^nx 
defined by 

$a(x) = $,(Ax) + In A' , 1 < i < n (45) 

are also solutions of eq. (p^ . This property is due to the scale behaviour of Newton's 
potential. 

In the spherically symmetric case the mean field equations (0) become ordinaries 
non-linear differential equations 



d^^i 2 d^i 



" 71- 

+ 4nnii e*^^^) = , 1 < i < n. (46) 



24 



Using the scale covariance of the mean field eqs.(0) by the transformation (^51), we can 
set 

^ * i'lt) MS) 

with new parameters 11^ = J'^ , 1 <i <n. In this way, the mean field eqs.fEBD become 
a reduced system of the form 

X^W + TX^W + mT. =0 ,l<i<n. (48) 

Let us find the boundary conditions for these equations. In order to have a regular solution 
at the origin we impose 

x'M = ,l<t<n. 
The system (H^) is invariant under the transformation 

A ^ A e", Xi Xi - 2a , 1 <i <n . 

Hence, we can choose 

Xi(0) = 

without loosing generality. As in the case of two kinds of particles, the remaining boundary 
conditions X2(0), Xn(0) are not independent from rji, 772, Vn- The normalization (|iT| ) 
of the densities of the n kinds of particles pi, p2, . . . ,Pn has to be imposed. We obtain 
from eqs.dD and (|^, 



- dx a;^ e^'(^) , 1 < i < n. (49) 
A Jo 



1 

Using the reduced mean field equations (^), it is straightforward to show from eq.(^ 
that 

mi ^ = -A xl(A) , l<i<n. 



Hence, 

x;(A) = ^x:(A) , \<%3<n. (50) 

From eq.(pOD we introduce new A-independent parameters 

vn ■ 

a = X.(A) - — Xi(A) , 2<2<n. 
mi 

These new parameters are only function of 77^, ■, rj^- Notice that the boundary 
conditions can be written as, 

X2(0)=C2, ... , Xn(0)=C„. 

The reduced mean field equations (^) become a single ordinary differential equation with 
its coefficients depending on the parameters C2, C„ , 

x'i'(A) + |x'i(A)+miX: — e^>^^(^)=0 (51) 
A ^ mi 

and the boundary conditions Xi(0) = 0, x'i(O) = 0. [Here, Ci = 0]. 
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Using eqs. (|43|) and (^Tf ) we can express the densities of the n kinds of particles in terms 
of the solution of eq.(pT|) 



Pi{R) 



R 



1 < i < n 



The thermodynamic functions are expressed as functions of the density of particles at the 
boundary /i, /a, ...,/„ depending on r/f , t]^, r/^. That is, 



f. - ^ ^ p;;rrXi(A) 



1 < i < n . 



We provide below the expressions for the free energy, the gravitational energy, the entropy 
and the equation of state. 



F-Fq 
T 

Ep 

S 
PV 



1=1 j= 

n 

3TJ2N^{f^-l] 



(52) 



i=l 



i=l 



(53) 



1=1 



where 



Fo = -J2NiT In 



i=l 



eV f mjT\ 2 



is the free energy and 



i=l 



V /m,T\^ 



V 27r 



is the entropy of the perfect gas with masses mi, m2,..., m„. Combining eqs.(^) and 
(0) yields the virial theorem 

3T 



T 



i=l 



5 Conclusions 



The self-gravitating gas with two kinds of particle has analogous qualitative properties 
to the self-gravitating gas with one kind of particles. Physical quantities like energy, free 
energy and entropy turn to be the sum of a term proportional to Ni plus another term 
proportional to N2 for large Ni, N2 and V provided Ni/V^ and N2/V3 are kept fixed. 
All physical quantities are expressed as functions of rji and 772. Instead of a critical line as 
for one kind of particles, we have here a critical line in the (r^i, 772) plane for the canonical 
ensemble . This line is associated to the Jeans instability. 

The equation of state exhibits a rim on this critical line [see fig.^. The thermodynamic 
functions exhibit a two-sheeted structure as functions of t]^ and 77^^. The branch points 
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are on the critical line. The specific heat is discontinuous and diverges there while the 
free energy is finite and continuous in the branch line. 

The local pressure and the local densities of particles are related by the same equation 
as in a perfect gas [see sec. 3.5]: 



P(r) 



NiT 



Pi(r) + 



V 



This can be explained by the dilute character of the self-gravitating gas in thermal equi- 
librium: N/V ~ N~'^ for N ^ oo. This dilution damps the effective interparticle 
interaction and allows a free particle behaviour. 

The particle distribution is inhomogeneous and scales with R with a Haussdorf dimen- 
sion d. The Haussdorf dimension d decreases for increasing r)i and r^l^. Its value on the 
critical line d = 1.6 .. . turns to be independent of the ratio Tj^/rj!^ implying an universal 
behaviour, d takes there the same value than for the canonical critical point with one 
kind of particles [ffl . 



.1 Appendix A 

The goal of this appendix is to compute the expression 



A=!^Jd'^ <l>i(x) e*^W + ^ /d'x $2(x) e*^W 



Using eq.(^ we obtain 

iVi 



m2N2 



rriiNi 



N2 , m2 

+ — (^2 CtlJ 

2 mi 



Using eq.(|TT]) in the spherical symmetry the first term of eq.(|5^) becomes 



Ni 
2 



d3x<|.i(x)(e*^W + ^^*^W 



miNi 



]\[-^ rRma. 



2r]i Jo 

Using eqs. (|T6|) and pO| ) we obtain 



_L aR^i{R)—{R 



X 



R-n 



dR^ dR'' 



Integrating by parts (|55D and using eqs. ([T6|) and (|56D we obtain 

m2N2 



J d^x<l>i(x) 



Ni 
2 



Ni , m2N, 



miNi 



Using eq.(p9|) the second term of (|5^ yields 



A^2 


/ "^2 \ 






02 ai 




2 


V mi J 





m2 
mi 



(54) 



(55) 



(56) 



(57) 



(5^ 
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Using eqs. (|57D and (|5^) we express A as 



We compute now 



/(A)= / dxx'[x[{x)f 



(59) 



Using eqs.(^) and (PT]), we derive the function 



and find 



(60) 



We integrate B' between and A. Integrating by parts the first term and using eg. (pO| 
we find 

The second term of eq. (|60|) yields 



-2/(A) . 

Using eq.(|19D the third term of (0) yields 

3A« + v^) . 

Hence, 

/(A) = -2A=^[e^^(^) + /i^e^^^^)] - Xir]^ + fi r]^f + 6A(r/f + /i^ t]^) 



Therefore, using eqs.(y) and (|28|) we obtain 



2\rj^ Jo 



dx X [Xi{x)] = Ni 



No 



R\ 



3(1 -/i)- 2(^1 +-"^2 
3(l-/2)-^(%f + r/,^' 



(61) 



Hence, the expression of A [using eqs. (|59D and (|6lD] is 

A = Ni 
+ N2 



^ln/i + 3(l-/i)-i(< + /xr/| 



2 

1 / 1 



-ln/2 + 3(l-/2)--(^^ V? + V^ 



.2 Appendix B 



We derive here the mean field equations from the hydrostatic equilibrium condition com- 
bined with the Poisson equation for a mixture of two ideal gases. 
The hydrostatic equilibrium condition 

VP(g) = -KiViPi(g) + m2iV2p2(g)] Vf/(g) , 

where P{q) stands for the pressure, combined with the equation of state for the ideal gas 
in local form 

P{q)=T[N,pM + N2P2m , 
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yields for the particle densities 

where pi and ^-^^g constants. Inserting this relation into the Poisson equation 

V^U{q)=4TTG [nil Ni pi{q) + 7712 N2 P2{q)] 
yields in dimensionless coordinates eq.(|np. 
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